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1 Introduction 



Gauge fixing is needed in manifestly gauge invariant lattice gauge theories for non- 
perturbative calculations of gauge dependent quantities such as fermionic and gauge 
field propagators. It became a useful tool for investigations of the vacuum structure 
in non-Abelian lattice gauge theories, too. 

However, one is faced with the non-uniqueness of the solutions of the commonly 
used gauge conditions, Lorentz, Coulomb gauge, etc. |]T|. Let us call the different 
solutions of the same differential gauge condition gauge copies. In practice, the 
existence of different gauge copies can strongly affect the numerical computation of 
gauge dependent quantities. 

For U{1) lattice theory the existence of gauge copies has been well established 
(see papers 0-0). It should be mentioned that the ambiguity also exists for the 
continuum U{1) gauge theory, if the latter is treated on a four-torus [0]. 

The existence of gauge copies raises the following questions. 

• Are there copies due to pure lattice artifacts? 

• What is the influence of gauge copies on physical quantities such as field 
propagators? 

• How should a unique gauge fixing prescription look like? 

It is the goal of this paper to provide answers for four-dimensional compact U{1) 
lattice gauge theory with periodic boundary conditions and within the Coulomb 
phase. 

Let us consider (lattice) gauge conditions of the type 

n n 

Y dl-^lm U,, = Y dj,-^ sin 9,, = 0, (1) 

where Ux^ = exp (iOxij.) and the lattice derivative dj^^ fx = fx — fx-p. ■ The case 
n = 4 represents the nonlinear (lattice) Lorentz gauge condition |^ - often called 
also Landau gauge -, whereas n = 3 leads to the Coulomb gauge. The differential 
gauge conditions are equivalent to finding stationary points of the functional F 

= ^E^-W ' = [ cos^.^ +cos^^,_^;^ ] (2) 

with respect to periodic gauge transformations (V denotes the number of lattice 
sites). Obviously, the requirement of a global maximization of expression (^ is 
stronger than the differential condition (|I]). 

The usual way to fix the gauge is as follows |^. The equilibrium configurations 
produced by Monte Carlo updating are independently gauge transformed in order 
to satisfy the gauge condition. In case of the gauge conditions (|I]) this is performed 
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by means of an iterative procedure which maximizes Fx{6) site by site. We 
shall call this algorithm standard gauge cooling. In order to improve convergence 
overrelaxation steps are usually built in. One can easily convince oneself numerically 
that standard gauge cooling for Lorentz and Coulomb gauge leads to various stable 
local maxima related to different F values. The following discussions will mainly 
apply to the Lorentz gauge. 

In a recent paper ||^, it was shown that in the Coulomb phase some of the gauge 
copies produce a photon propagator with a decay behaviour inconsistent with the 
expected zero mass behaviour. The origin of these 'pathological' gauge copies was 
not clarified, but a recipe to avoid this problem was proposed. According to this, 
one has to generate several gauge copies by applying standard gauge cooling to 
randomly gauge transformed equilibrium configurations. Then choosing the gauge 
copies with maximal F{6) values significantly improves the photon mass estimate 

i- 

In the following we want to go a step further. We are going to show that the 
'bad' gauge copies are related to the existence of pairs of Dirac sheets, and we 
propose an algorithm which removes them in a straightforward way. 



2 Lattice actions 

We consider U{1) lattice gauge theory on the torus . In our calculations, we 
first apply the standard Wilson action (WA) 

swm = p-T.i^-^o^^p) ' (3) 

p 

where plaquette angles 9p = 9^^^^ are given by 6*^. = 9^^ + 9^+f,^ ^ - 6*^.+^, ^-9^^ . 
The plaquette angle 9p can be split up 

9p = 9p + 27mp, -n<9p<7i, np = 0,±l,±2, (4) 

where the gauge invariant part of the plaquette angle, 9p , describes the electro- 
magnetic flux through the plaquette, and np is the number of Dirac strings passing 
through it. Let us call plaquettes with np ^ Dirac plaquettes (DP). 

Monopole charges Kc in elementary 3D-cubes c are then defined by the net 
number of strings entering or exiting these cubes. The worldlines of these Abelian 
monopoles on the dual lattice close either within the lattice volume or by periodic 
boundary conditions. The dual integer valued plaquettes 

m,,,, = -£^.p.n,,p. (5) 

form Dirac sheets bounded by the worldlines of monopole-antimonopole pairs or 
closed within the lattice volume or due to periodic boundary conditions. 
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In a recent paper |T0[, we showed that the monopoles (as well as negative pla- 
quette values) have to be understood as lattice artifacts. Their suppression removes 
the phase transition of the standard Wilson theory at /3 ~ 1 and as a consequence 
provides a uniform Coulomb phase in the whole positive (3 range (see also ||11|| ). 
Therefore, we will discuss also the modified compact action (MA): 

Sma = Sw + >^k -^2 \ +Az ■ H ^ sign(cos 6'p)) , Xk,Xz = oo . (6) 

c p 

Note that this modification does not influence the classical continuum limit and 
does not change any perturbative aspects. It may serve as an alternative for non- 
perturbative QED investigations carried out so far with the non-compact gauge 
action. 



3 Gauge fixing and Dirac sheets 

Most of our calculations were done on 6^ and 12 ■ 6^ lattices, for some exploratory 
calculations also on 20^ and 28 ■ 14^ lattices, respectively. 

First we applied ten random gauge transformations to each Monte Carlo equi- 
librium configuration. On each of these versions, we carried out standard gauge 
cooling with the stopping criterion 

Max I J2 dl'^ sin e,^\ < lO'' and ^ E I E ^i"^ sin e^^\ < 10"^ . (7) 

We determined the average number (Nim) of different local maxima Fmax{d) of 
the functional (^ found for each equilibrium configuration. Numerically, two local 
maxima are understood to be identical if the difference of corresponding Fmaxi^) 
values is less than e = 10~^ . (To make sure that this recipe is well-defined, we 
varied the value of e and found that for e = 10~^ and 10~^ , the change in {Nim) 
was less than the statistical errors.) For the Wilson action in the confinement phase, 
we found {Ni^) ~ 10 ; i. e. practically all gauge copies produced are related to 
distinct values F^axi^^) - At /3 ~ 1 , this average drops sharply down to {Ni^) — 3 
and approximately does not depend on (3 within the Coulomb phase. For the 
modified action, the value of (Nim) is seen to be independent of the coupling in 
the whole interval of /3's considered and coincides with that for the Wilson action 
in the Coulomb phase. 

The functional dependence of {Nim) = f{P) for both actions turned out 
to look very similar to the behaviour of the average number of Dirac plaquettes 
{Ndp) = giP) after standard gauge cooling. Therefore, we guess that both the 
values of Fmax{d) and the number of different local maxima are correlated with 
{Ndp) ■ 

Another observation made in both the Coulomb and the confinement phase (for 
Wilson action) proves this idea. There is a strong correlation of the local quantity 
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Fx{d) with the number of Dirac plaquettes among the 24 plaquettes sharing the 
site X . 

By a more detailed investigation of the number and the location of Dirac pla- 
quettes in the Coulomb phase, we observed pairs of Dirac sheets closing by periodic 
boundary conditions. In the extreme case of a minimal surface this is a pair of 
periodic Dirac sheets occupying parallel planes and having corresponding Dirac 
plaquettes with opposite signs. In general the surfaces fluctuate around this ex- 
treme case. The minimal number of Dirac plaquettes in such a lattice gauge field is 
equal to 2 ■ ■ L^, , where L^j^.Ly , fj,,^ = 1, 4 denote the corresponding linear 
extensions of the lattice. 

In contrast to monopoles, gauge transformations can move and deform Dirac 
sheets. Therefore, we should be able to remove double periodic Dirac sheets - called 
below simply double Dirac sheets (DDS) - by gauge transformations. However, it is 
worth noting that this property is in contrast to that of single periodic Dirac sheets 



as investigated e.g. in |12 



In the following we shall discuss gauge fixing with a special focus on the role 
of DDS. To each equilibrium configuration obtained originally from a cold start, 
we apply the standard gauge cooling procedure only once. Fig. la shows the 
resulting distribution P{Nj:,p) for the modified action on a 12 ■ 6^ lattice for 
f] = 1.1 . There are only lattice configurations either with N^p = or 
Ndp > 2 ■ 6^ = 72 , the latter attributed to DDS. For the standard Wilson action 
at (3 = 1.1 the corresponding distribution looks less pronounced because of 
monopoles, which are still allowed in the Coulomb phase. But in this case, another 
plot makes the situation transparent again. On each configuration obtained after 
gauge fixing, we measure the number of Dirac plaquettes separately for every two- 
dimensional plane orientation and select the one with the maximal number N'^p^ . 
Fig. lb shows the corresponding distribution for the N^p^ values. There are 
two nice distinct peaks in the distribution - related to configurations either with or 
without DDS. 

There are no DDS (or even single ones) in the ^ — //-plane provided N^p^ < 
■ is satisfied. We can use this criterion for dividing the gauge copies into 
two distinct classes: those which definitely do not contain DDS and those where 
DDS cannot be excluded. Fig. 2 shows a scatter plot for the local maximum values 
Fmax{G) found for 1000 equilibrium configurations after standard gauge cooling 
for the Wilson action (WA) at (3 = 1.1 . One can see that the crosses (copies 
without DDS) lie systematically higher than the full squares (copies presumably 
with DDS). The corresponding average values are indicated by broken lines. For 
the modified action (MA) at the same (3 , the analogous plot looks similar, but 
with the averages shifted to larger Fmax values. 
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4 The photon correlator 



Let us now consider the (lattice) photon propagator. The corresponding correlation 
functions r(r;p) are defined by 

TiT;p) = {^*it + T;p)-^t;p)) , <|.(t; p) = ^ 6""^^ x) , (8) 

X 

with 

• the gauge dependent link operator 

=sin^,, , 1 = 1,2,3 (9) 

or 

• the gauge invariant plaquette operator for the state 

= sin^,,,„ t,j = 1,2,3 (10) 

For no n- vanishing momentum the latter state has an overlap with the 1 
photon state |TB . 



For the gauge invariant operator, the correlator behaves as 

r(r;p-) 
r(0;p) 



A- 



exp(-r-E(p))+exp(-(L4-r)-E(p)) ]+... , 0<A<1 (11) 



where dots correspond to higher state contributions. In the massless case, the 
energy E{p) is given by the lattice dispersion relation: 

sinh2- = f]sin2^ . (12) 
2 i=i 2 

Provided the Lorentz gauge is properly fixed, the transverse photon propagator in 
the Coulomb phase is expected to have the same simple form eq. ([TT]). We used the 
momenta p = ^ ■ {0; k; 0) , k = 1,2 in computing the correlators and accordingly 
the gauge dependent operator (y) with i = 1,3 as well as the gauge invariant one 
(0)with (z,j) = (l,2),(2,3). 

Our results for all the correlators are plotted in Fig. 3a in case of the Wilson 
action at jS = 1.1, i. e. in the Coulomb phase. One can see that the propagator 
obtained from operator (^ and averaged over all configurations without taking 
care of the presence of DDS (crosses) clearly deviates from the theoretical curve, 
eqs. (|TI| , 1^). Full circles correspond to an average over such gauge copies which 
presumably possess DDS. The deviation from the theoretical curve becomes even 
more pronounced in this case. But the correlator obtained after averaging only 
over those gauge copies which definitely do not contain DDS nicely agrees with 
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the theoretical expectation. We show also the correlator for the gauge invariant 
operator eq. ( |T0|) (squares). For the Wilson action its points lie somewhat below 
the theoretical curve. This is a known effect |Ty] which is due to the fact that the 
overlap coefficient A for the plaquette operator is less than one, i.e. there is a 
non-negligible admixture of higher states. 

Fig. 3b shows the analogous results for the modified action at f3 = 0.9 . Here 
we obtain a perfect agreement between the photon correlator measured with both 
the operators (H) - if configurations with DDS are excluded - and (|T0D and the 
theoretically expected curve. 

An analogous investigation for the transverse photon correlator in the Coulomb 
gauge produced similar results. 

Thus, we conclude that the presence of DDS as gauge and lattice artifacts spoils 
the correct determination of the photon mass from the gauge dependent operator 



5 Improved gauge fixing algorithm 

A Polyakov loop investigation will be helpful to propose a gauge fixing algorithm 
allowing to remove DDS. 

For a given configuration on a lattice with linear extensions , fi = 1, ••■,4, 
we define the average Polyakov loop along direction fj. 

^=^E^^-M = ra■e^'^^ -^<^,<^, (13) 

where = V/L^ . In Fig. 4 for the Wilson action at /5 = 1.1, we show the 
scatter plot of values in the complex plane for 400 Monte Carlo equilibrium 
configurations investigated after standard gauge cooling. The distribution has the 
form of a circle, as it should be. But the points corresponding to DDS orthogonal 
to the //-direction are not distributed with rotational symmetry. They occupy only 
the left part of the circle, i.e. Re(P^) < . 

It is worth noting that gauge transformations which are implemented in the 
standard gauge cooling procedure are periodic by construction. Each local 
value cannot change under standard gauge cooling transformations, and therefore, 
their distribution is also invariant. 

This suggests how to avoid gauge copies with DDS. We carry out nonperiodic 
gauge transformations 

Ua:,,^ 9*-UxfJ- 9x+f, = Uxfj-e''^, ^^ = 6^""^'^ (14) 

or equivalently 
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which turn all P^, fi =1,...,4 into real positive numbers before applying standard 
gauge cooling. Therefore, we put 

c^ = -i^^L/L^„ /i=l,...,4. (15) 

As is evident from Fig. 4, there are no gauge copies with DDS among configurations 
with large positive real part of P„ . This gives us confidence that after a gauge 



transformation of the type (0), which we will call preconditioning step, subsequent 
gauge cooling should avoid local maxima with DDS. Our computations (with around 
30000 measurements for the same f3 and action as in Fig. 4) confirmed this 
conjecture. No gauge copies with DDS were found. 

The reader will be convinced also by Fig. 5. Starting for every equilibrium 
configuration with the preconditioning step, we look after each subsequent gauge 
cooling sweep for the layer with the maximal number of Dirac plaquettes N^p^ as 
explained above. N^p^ has been plotted for all 400 configurations investigated as a 
function of the number of gauge cooling sweeps Ngg^ . On a 12-6^ lattice, already 
after 12 sweeps no configurations with N^p^ > 36 were left. We checked the 
algorithm described above also with larger lattices ( 14'^ ■ 28 and 20^ , respectively) 
in the Coulomb phase. There were no DDS. 

It is important to stress that the above procedure avoiding DDS already allows 
us to reproduce the correct photon propagator despite the fact that it does not 
yet provide unique gauge fixing. The global maximum is not really necessary for 
solving the photon mass problem. 

However, it appears that a combination of nonperiodic gauge transformations 



(p^) and gauge cooling permits to find the global maximum of F. The recipe looks 
as follows. After the preconditioning step, we carry out gauge cooling sweeps until 
DDS are removed (around 10 maximization and overrelaxation sweeps for the lattice 
sizes discussed). In a third step, we maximize F{9) with respect to nonperiodic 
gauge transformations (|T1) . Thus, the global phase shifts ,/i = 1, ...,4 are fixed 



such that 

F,{0)^^T.U.,^\F,\-e^^-^\F,\, /X =1,...,4 (16) 

become real and positive. Since this maximization generally violates the local gauge 
condition (|l|) we mix it with further gauge cooling iterations. This procedure can be 
repeated until the changes in F{6) during the nonperiodic gauge transformations 
vanish within machine precision. We checked the whole recipe by applying it again 
to 10 randomly gauge transformed versions of each Monte Carlo equilibrium field. 
Only for 0.5% of all equilibrium configurations we still found gauge copies. They 
are, however, related to each other by gauge transformations Qx significantly 
differing from unity only in a very narrow neighbourhood of a single site! 

The algorithm proposed here applies also to the Coulomb gauge in an analogous 
manner. 
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6 Conclusions 



Let us summarize our findings. 

• We have convinced ourselves that gauge copies which possess double Dirac 
sheets (DDS) give rise to a wrong behaviour of the transverse photon propa- 
gator in the Coulomb phase. 

• The transverse photon propagator obtained by averaging over gauge copies 
without DDS shows an agreement with the expected theoretical curve which 
for the standard Wilson action near to the deconfinement phase transition is 
even better than for the gauge invariant plaquette correlator. 

• We proposed and checked an algorithm for Lorentz (and Coulomb) gauge 
fixing which avoids gauge copies with DDS. 

• Nonperiodic gauge transformations were shown to allow unique gauge fixing 
to the global extremum of the gauge functional F{6) in almost all cases. 
This might be of importance for fermion propagator investigations now under 
way. 

We did not consider the confinement phase in case of the standard Wilson action. 
Because of the existence of monopole excitations, a unique gauge fixing is more 
involved than for the Coulomb phase. Since this is even more important for non- 
Abelian gauge theories, we hope to come back to this point very soon. 
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Note added in proof 

After having written up this letter, we became aware of the recent paper by M.I. Po- 
likarpov. Ken Yee and M.A. Zubkov, preprint LSU-431-93, hep-lat/9305019 , where 
some of the conclusions about the role of Dirac strings in lattice QED are drawn, 
too. 
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Figure captions 



Fig.l. Distributions of the number of Dirac plaquettes PiNop) on a 12-63 lattice 
at (3 = 1.1 for modified action (MA) (a) and for Wilson action (WA) (b). In the 
latter case, only those plane orientations are selected which contain the maximal 
number of Dirac plaquettes NJ^p^ . 

Fig. 2. Fmax{(^) for different configurations for WA on a 6^ lattice for (3 = 1.1 . 
Crosses (squares) correspond to configurations without (presumably with) double 
Dirac sheets. Broken lines are the corresponding averages for the two classes of 
configurations. 

Fig. 3. Photon correlators T{t;p) / T{0;p) with momentum p = ^ ■ (0; 1; 0) on 
12 ■ 6^ lattices at /3 = 1.1 for WA (a) and at /3 = 0.9 for MA (b). Open (full) 
circles correspond to transverse correlators calculated on the configurations without 
(with) double Dirac sheets. Crosses show the transverse correlators averaged over 
all configurations. Squares correspond to plaquette-plaquette correlators. The solid 
line represents the expected zero-mass behaviour (|1^). 

Fig. 4. Scatter plot for Polyakov loop values P^, /i = 1, 2, 3 in the complex 
plane after standard gauge cooling (WA, /3 = 1.1 , 12 ■ 6^ , 400 measurements). 
Squares (crosses) correspond to the existence (non-existence) of double Dirac sheets 
in planes orthogonal to /i-direction. 

Fig. 5. N^p^ versus number of gauge cooling sweeps Ngsw after preconditioning 
(WA, f3 =1.1 , 12 ■ 6^ ,400 measurements). 
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